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Abstract 

Nonstandard analysis is an area of modern mathematics which 
studies abstract number systems containing both infinitesimal and 
infinite numbers. This article applies nonstandard analysis to de- 
rive jump conditions for one- dimensional, converging shock waves in a 
compressible, inviscid, perfect gas. It is assumed that the shock thick- 
ness occurs on an infinitesimal interval and the jump functions in the 
thermodynamic and fluid dynamic parameters occur smoothly across 
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this interval. Predistributions of the Heaviside function and the Dirac 


delta measure are introduced to model the flow parameters across a 
shock wave. The equations of motion expressed in nonconservative 
form are then applied to derive unambiguous relationships between 
the jump functions for the flow parameters. 


1 Introduction 

Currently, no known mathematical theory provides a general framework to 
study all the partial differential equations of mathematical physics. The lack 
of such a comprehensive theory means that each type of physical problem, 
and hence, each type of equation, requires a different analytical or numerical 
technique for its solution. Functional analysis has been applied to develop 
powerful techniques for analyzing and constructing generalized solutions to 
families of linear partial differential equations. However, the methods of 
functional analysis are limited in studying nonlinear partial differential equa- 
tions, because such techniques are based on linear mathematical structures. 
Since linear mathematical structures do not describe the simplest nonlin- 
earity (multiplication), it is not surprising that functional analysis has not 
yielded a comprehensive framework for studying generalized solutions of par- 
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tial differential equations. 

Many approaches have been used to attempt to define a useful multiplica- 
tion operation on spaces of functionals applied in linear analysis to study dif- 
ferential equations. The motivation for multiplying functionals comes from a 
plethora of applications that occur in modeling physical problems: construc- 
tion of generalized solutions of certain partial differential equations requires 
the multiplication of linear functionals. However, a fundamental theorem of 
Schwartz [1] shows that the standard algebraic product cannot be defined 
consistently in spaces of functionals, such as the distributions or generalized 
functions; all such methods will produce ambiguous results leading to log- 
ical and arithmetical contradictions. Because standard multiplication does 
not produce meaningful results, any product of generalized functions defined 
must be weaker than standard multiplication. 

Over the last thirty years, E.E. Rosinger [2 - 6] and J.F. Colombeau [7, 
8], working independently, have developed a theoretical framework for quo- 
tient algebras of generalized functions to analyze nonlinear partial differential 
equations. Their approach extends classical linear mathematics to nonlinear 
mathematics by defining multiplication on sets of generalized functions by 
embedding these linear spaces into abstract algebras of nonlinear generalized 
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functions with weak operations of multiplication. The potential applica- 
tion of a theory of nonlinear generalized functions enabling the analysis and 
approximation of solutions to nonlinear partial differential equations is enor- 
mous; such a theory could impact almost all areas of science and technology. 

Although it is probable that the existing theories of quotient algebras of 
generalized functions will yield important insights into the analysis and com- 
putation of certain classes of nonlinear held equations, significant limitations 
of these theories are known. Because any logically consistent operation of 
multiplication of generalized functions must be much weaker than the clas- 
sical multiplication of numbers and functions, such a product of generalized 
functions will always yield large collections of abstract objects representing 
weak solutions of the problem under investigation. Basically problems with 
no classical weak solutions, such as the product, 5 ■ H, are converted into 
problems with uncountably many generalized solutions. The fundamental 
question then becomes how to determine a physically observed weak solu- 
tion from an infinitude of possible weak solutions. 

To approximate solutions, Colombeau (see References [7, 8]) constructed 
an explicit theory of nonlinear generalized functions. This theory has been 
applied to analyze and compute solutions to nonlinear partial differential 
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equations. The logic of solving a differential equation in a space of Colombeau 
generalized functions follows the standard approach of functional analysis: 
a method is established for generating sequences to approximate solutions 
within a collection of nonlinear generalized functions and then a compactness- 
type argument is applied to show that the approximating sequences actually 
converge to limits within the space of approximating objects. A resulting 
limit then represents a solution to the equation. A key difficulty in applying 
a space of Colombeau generalized functions is determining a specific represen- 
tation for the products of generalized functions occurring in the differential 
equation. As Colombeau [9] and his coworkers [10, 11] have shown, differ- 
ent physical problems require different distinct representations for a product 
of generalized functions such as S ■ H . The representation of the product is 
found by adding more information about the problem. Roughly speaking, the 
solution of a concrete problem via nonlinear generalized functions requires: 
the equations modeling the physics, the initial-boundary values, and the ad- 
ditional problem-specific physical information required to fix the product of 
generalized functions. For example, in continuum shock wave problems, the 
additional information is obtained by requiring the product of generalized 
functions to satisfy the standard shock wave jump conditions. 
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A number of authors, including, Oberguggenberger and Todorov [12], and 
Hoskins and Sousa Pinto [13], have shown that versions of the Colombeau al- 
gebras of nonlinear generalized functions may be constructed using nonstan- 
dard analysis. Nonstandard analysis is a relatively new area of mathematics, 
which emerged from basic research in mathematical logic. The subject was 
discovered in the early 1960s by A. Robinson [14], The main contribution 
of nonstandard analysis to mathematics is the extension of the real numbers 
R to the hyperreal numbers *R that contain infinitely small (infinitesimals) 
and infinitely large numbers. The infinitesimals in a hyperreal number sys- 
tem have the algebraic properties of standard numbers, which justifies the 
formal manipulations of infinitesimals that engineers and physicists often 
use. Moreover, the existence of distinct, infinitely large numbers that can 
be manipulated and used like finite numbers to solve problems provides a 
powerful new tool for applications. Nonstandard methods may be used to 
construct algebras of nonlinear generalized functions, or nonstandard analy- 
sis may be used directly to study specific problems involving multiplication 
of generalized functions. The second approach is considered here. 

In this article, jump conditions are derived using nonstandard analysis for 
one-dimensional converging shock waves in a compressible, inviscid, perfect 
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gas. A class of nonstandard Heaviside functions is applied to model the 
microstructure of the shock wave jump. The term microstructure refers to the 
distinct functions, which represent the flow variables such as specific volume, 
pressure, and velocity across a shock wave. A shock wave is governed by the 
standard conservation equations written in primitive form. A shock jump 
is then modeled to occur over an infinitesimal interval with nonstandard 
jump functions. The nonstandard jump functions applied in the subsequent 
analysis are defined as follows: for a standard jump function 4>(x) with jump 
[0] = (p r — <t>i at x = 0, define a nonstandard jump function by 

*(p(x) = 0/ + [<f>]*H(x), (1) 

where *H(x ) is a piecewise differentiable, nonstandard Heaviside function. 
Following Colombeau [9], each flow parameter is assumed to have a distinct 
nonstandard jump function associated with it at the shock. This means that 
each flow variable may have a different nonstandard Heaviside function or 
microstructure across a shock wave. 

ft is shown that if the various Heaviside functions for the flow param- 
eters have their jumps located on the same infinitesimal interval, then the 
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governing equations expressed in nonconservative form yield the relations 
between these Heaviside functions unambiguously. Strictly speaking, the 
mathematics of the product of nonstandard generalized functions is unam- 
biguous; therefore, the use of nonstandard Heaviside functions removes some 
of the mathematical difficulties associated with products such as 6 ■ H . 

The present work was motivated by the authors’ interest (see Baty et ah, 
[15, 16]) in understanding how to add physical information to specify prod- 
ucts of generalized functions contained in differential equations with model 
problems with well-defined regions of solution discontinuity. Nonstandard 
analysis is applied to build explicit examples of the infinitesimal microstruc- 
ture for converging shock waves. If the existence of the nonstandard Heav- 
iside functions is assumed, the present development greatly simplifies the 
analysis and manipulation of products of generalized functions. 

2 Nonstandard Generalized Functions 

Nonstandard analysis studies the extension of number systems and function 
spaces to quotient spaces, which contain idealized elements that are infinitely 
large and infinitely small. For example, the real numbers, R, may be ex- 



tended to a hyperreal number system *R that is defined as a quotient space 
resulting from applying an equivalence relation ~ on the set of sequences of 
real numbers, R s , indexed by the natural numbers N 

*R = Rs/ ~ . (2) 

The equivalence relation ~ is defined by selecting an ultrafilter on N. A 
resulting nonstandard set *R may be shown to be a linearly ordered field. 
The role of the equivalence relation in Equation (2) is to add the idealized 
elements. 

For the present analysis, let S denote the real numbers or a set of functions 
such as the space of locally integrable functions, Li oc (R). Let S' R+ be the 
set of nets (/ e ) in S with parameter e G R+. The nonstandard extension 
of S will then be defined to be the quotient space *S = S R+ / ~ for a 
fixed equivalence relation. The use of nets provides a generalization of the 
definition of Equation (2). Further details of nonstandard extensions are 
available in many textbooks, articles and proceedings; for example see Hurd 
and Loeb [17] and Arkeryd, et al. [18]. 
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2.1 Predistributions of Generalized Functions 


In the early development of nonstandard analysis, hyperreal techniques were 
applied to generalized functions, Robinson [14], with application to quantum 
theory, Kelemen and Robinson [19]. Robinson noticed that any regularizing 
sequence of smooth functions used to represent the Dirac delta measure could 
be replaced by a single function with hyperreal parameters. A nonstandard 
Dirac delta measure results for any standard function g that satisfies 



by defining the internal function, 


*5(x ) = Dg(fhr), 


( 4 ) 


where 12 is any infinite hyperreal number. To see how the internal function 
(4) acts like the delta measure, perform the integration 
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where e = e(0) is an infinitesimal constant depending on If, 0 G D( R) is 
an arbitrary standard test function, and the equivalence relation « means 
infinitesimally close. The result follows by taking the standard part of Equa- 
tion (5). 

The power of the linear theory of generalized functions comes from the 
fact that many distinct regularizing sequences of functions may be applied to 
produce unique functionals, such as the delta measure. It is well known that 
the standard theory of generalized functions may be extended to a nonstan- 
dard theory of generalized functions, Laugwitz [20], Richter [21] and Todorov 
and Vernaeve [22] . In standard constructions, the delta measure results from 
abstract limits of regularizing sequences; on the other hand, in nonstandard 
constructions, the functional whose standard part produces the same result 
as the delta measure can be represented by uncountably many internal hy- 
perreal functions with distinct microstructures. These internal functions are 
called predistributions and are defined here for the space of locally integrablc 
functions *Li oc (R). A function / G *Li oc (R) is a called a predistribution if 
for each n G N 0 it satisfies the following: i) /-T is piecewise differentiable, ii) 

^ tE L loc (R), and iii) there exists a generalized function T G D'{ R) such 
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that 


/ ( ”>(x)>(x ) Ax ss (-l) n (T,<t> M ), 


(6) 


for all test functions (j) € D( R) and all n E Nq. 

The analysis of shock wave jump conditions also requires the use of non- 
standard Heaviside functions. Like the Dirac delta measure, the standard 
Heaviside function can be represented by uncountably many functions with 
distinct microstructures (predistributions of the Heaviside function). To in- 
troduce a predistribution of the Heaviside function, let and 62 be two 
infinitesimal numbers, at least one of which is not zero, ~ 62 and €\ < €2- 
An internal Heaviside function *H(x) is defined as follows: 


0 if x < €1 


*H(x) 


*h(x) if ei < x < €2 ' 


1 if x > e 2 


(7) 


In Equation (7) *h(x) is assumed to be a piecewise differentiable function 
contained in *Li oc (R), not necessarily monotonically increasing, which sat- 
isfies 

*h(e 1) « 0 and *h(e 2 ) ~ 1. (8) 
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As an example, the following function can be used for *h(x) when x is in the 
open infinitesimal interval (0, e): 


*h(x) = 1 — exp(— fix), (9) 

where is any infinite hyperreal and e = e(0) is any infinitesimal hyperreal 
such that fie is an infinite hyperreal. Combining Equations (7) and (9) yields 
the nonstandard Heaviside function 

if x < 0 

exp(— Qx) if 0 < x < e ■ 
if x > e 

Notice that *H(x) is a differentiable function on (0, e) and that the deriva- 
tive is the internal function 

' 

0 if x < 0 

H (x) = < y2exp(— fix) if 0 < x < e > (■*■■*■) 

0 if x > e 
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which is a predistribution for the Dirac delta. For all standard test functions 


4> 6 D{ R) it may be shown that 



*H\x)*<t>{x) dx 


m, 


( 12 ) 


following Equation (5). For more detail about nonstandard Heaviside func- 
tions and their properties see Baty et al. [16]. 


2.2 Products of Generalized Functions 

The Schwartz theorem showing that generalized functions cannot be multi- 
plied is a manifestation of a basic structural limitation of algebraic mathe- 
matics. Following Rosinger [5], consider a mathematical structure (J-'.-.'D) 
consisting of a collection T of functions defined on R that may exhibit a 
finite number of localized discontinuities. The standard Heaviside function 


H{x) = { 


0 if x < 0 


1 if x > 0 


(13) 


is contained in T . Let • and D denote the standard multiplication and 
differentiation operators acting on the set T . 
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From Equation (13) it follows that 


H m (x) = (H ■ H ■ H ■ . . . ■ H)(x) = H{x), (14) 

for all m > 1. Then applying the derivative operator in the sense of general- 
ized functions to Equation (14) yields 

mH m -\x)DH(x) = DH(x), (15) 

which implies that 

mH rn ~ 1 (x) = 1 and so H m ^ 1 (x) — — , (16) 

m 

for all m > 1. But Equation (16) contradicts the definition of the standard 
Heaviside function (13). This contradiction occurs because the Heaviside 
function has a jump at x — 0. Rosinger has shown that the mathemati- 
cal structure (T.-.V) is over-specified and discontinuity, multiplication and 
differentiation cannot be combined in the steps of a single derivation. The 
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over-specification of (JF, - ,T>) may be repaired by requiring that 


H 2 (x) = ( H-H)(x ) ^ H(x). 


( 17 ) 


From the definition of the nonstandard Heaviside function, Equation (7), 
it follows that 


* rj2 


0 


if x < E\ 


H ( x ) \ *h 2 (x ) if ei < x < e 2 ’ 


(18) 


if x > e 2 


which implies that 


*H 2 (x) ± *H{x) 


(19) 


because *h 2 (x) ^ *h(x) for < x < e 2 . That is, giving the Heaviside 
function *H(x) an inhnitesimal microstructure removes the over-specification 
which leads to the contradiction of Equation (16). Equation (19) generalizes 
for all differentiable nonstandard Heaviside functions to the result that 


*H m (x) ^ *H{x) 


( 20 ) 


for all m > 1, removing the contradictions obtained by Equation (14). Also 
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notice, that differentiating Equation (18) yields 


0 if x < €1 


*H(x)*H'(x) = < 


*h(x)*h'(x) 


if e 1 < x < e 2 ' 


( 21 ) 


0 if x > e 2 


As an example of a product of generalized functions, consider <5 ■ H. The 
multiplication of the delta measure with the Heaviside function occurs in the 
analysis and computation of shock waves if the governing equations of motion 
are written in nonconservative form; 5 ■ H comes from the nonlinear terms, 
such as the convective acceleration term uu x , for a discontinuity in the held 
variable, u. In the standard theory of generalized functions, the product 5-H 
does not yield a well-defined functional, because the results depend on the 
representations of 5 and H. To see this behavior, choose the predistributions 
of Equation (10) and 

*S(x) = — • sech(fhr), (22) 

7 r 

for the Heaviside function and Dirac delta respectively. Now, consider the 
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functional acting on any standard test function 0 e D( R) 


n 


*8(x)*H(x) (j)(x)dx = / — sech(flr)(l — exp(— flx))*0(x)dx. (23) 

Jo 71 


Integrating Equation (23) then gives 


r°° Q f i in 2 

/ — sech(fhr)(l — exp(— flx))*<f>(x)dx ~ 0(0) < 

/ 0 7T (2 7T 


(24) 


If other predistributions are selected for 6 and H, the numerical value mul- 
tiplying 0(0) in Equation (24) will, in general, be different. 

In the numerical approximation of shock problems, the use of noncon- 
servative forms of the equations of motion produce shock wave speeds that 
depend on the numerical scheme and not the physics of shock propagation. 
Nonconservative numerical simulations of one-dimensional shock problems 
are analogous to the result of Equation (24): the specific numerical values 
are fixed by the choice of the approximations for 5 and H . 

Next, consider the product 6 ■ H, but now assume 5 — H' . Then for any 
predistribution representing the delta measure, H' defined on the interval 
(— e, e), choose hyperreal numbers fl, e = e(f2) so that is infinite, e is 

infinitesimal and fie is an infinite hyperreal. Then for any standard test 
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function 0 G _D( R), the nonstandard functional equation for <5 ■ H implies 


' *5(x)*H(x)*(j)(x)dx = / * H' (x)* H (x)* <fr(x)dx, 


( 25 ) 


and 


*5(x)*H(x)*4>(x) da: « -0(0), 


(26) 


independent of the order of the operations used to evaluate the integral. 

Equation (26) is consistent with the results obtained by the standard 
theory of generalized functions. This special case of the product <5 • H may 
be computed using only the function H and multiplication and is not over- 
specified, since 5 ■ H is not differentiated. The differentiation of the Heaviside 
function is avoided by applying integration by parts. To this end, consider 
the standard version of Equation (25) 



S(x)H(x)4>(x) dx 




4>(x)dx, 


(27) 


which assumes that H{x) = H 2 (x). Equation (27) then yields 



-H(x) ) <j)(x)dx 


1 

2 



— 0(x ) dx 
ax 



(28) 
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producing the same result as the standard part of Equation (26). This exam- 
ple is analogous to integrating the governing equations of motion modeling 
shock wave propagation in conservative form: the predicted shock wave speed 
will be independent of the numerical approximation. Notice that the stan- 
dard calculation performed in Equations (27) and (28) depends on the order 
of the algebraic and analytical steps. Changing the order of the operations 
can produce an over-specified problem with ambiguous results. On the other 
hand, the nonstandard result of Equation (26) is independent of the order of 
operations. 

3 Shock Wave Jump Conditions 

3.1 Equations of Motion 

The equations governing the motion of a one- dimensional, compressible, in- 
viscid gas are given by 

t\jU 

(In p) t + -u(ln p) r + u r + — = 0, (29) 

u t + uu r + up r = 0, (30) 
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( 31 ) 


, I l/Ul V 

e t + ue r + up{u r H ) = 0, 

r 

where zz is the specific volume, p is the density (= ^),u is the velocity, p is the 
pressure, and e is the specific energy. The notation for partial differentiation 
is simplified to the subscript notation £ r = |£. Here k — 0,1,2 represents 
Cartesian, cylindrical and spherical coordinate systems, respectively. The 
conservation of mass, momentum and energy Equations (29), (30), and (31), 
represent three equations in terms of four physical variables. To close the 
system of equations two equations of state must be added. The equations 
of state relate the thermodynamic variables p, p, e, and T. Assuming a 
calorically perfect gas, the energy equation may be shown to be 

— — A(pv)t + u(pu) r \ + L>p(u r + — ) = 0, (32) 

7 — 1 r 

where 7 is the ratio of specific heats. Developments of Equations (29), (30), 
and (32) and discussions of the equations of state may be found in Zel’dovich 
and Raizer [23] and Thompson [24], 
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3.2 Jump Conditions for Normal Shock Waves 

To motivate the physical assumptions used to derive jump conditions for 
converging shock waves, consider the motion of a normal shock wave (k = 0) 
propagating through a perfect gas. On either side of the shock wave the 
specific volume u, velocity u, and pressure p of the flow are assumed to be 
constant. Since the values of the field variables are assumed constant on 
either side of the shock wave, the shock wave does not accelerate and the 
shock speed c is constant. Therefore, the characteristics for the shock wave 
are straight lines in the space-time domain. Along the characteristic lines 
the governing equations reduce to ordinary differential equations, so the field 
variables, u, u, and p , may be assumed to have the following form across the 
shock wave 


KO = vi + 

( 33 ) 

«(0 =ui + [u]K(£), 

( 34 ) 

P(0 - Pi- \p\L(£)f 

( 35 ) 


where £ = x + ct is a characteristic line, c is the shock speed, and where 
[0] = 0 r — 0/, with the subscripts r and l referring to the right (downstream) 
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and left (upstream) conditions across the shock, respectively. The functions 
H, K, and L are assumed to be predistributions of the Heaviside function 
as defined in Section 2.1. (Here and in the following development the *- 
notation is suppressed.) For an arbitrary fixed infinitesimal e, each one of 
these Heaviside functions is taken to have its jump contained on the same 
arbitrary infinitesimal interval (0,e). 

Now, substituting the jump functions (33), (34), and (35) into the equa- 
tions of motion (29), (30), and (32) yields 


u[v]H' — v[u\K' = 0, 

u[u]K' + v\p\L' = 0, 
au(p[u\H' + v\p\L') + vp[u\K' = 0, 


(36) 


(37) 

(38) 


where 


u = (ui + c) + [u]K and a 


1 

7-1' 


(39) 


Recasting Equations (36), (37), and (38) as a first-order system of ODEs 
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generates the formal equation 


- 


~ 

t 

- 

[u]u — [u]u 0 


H 


0 

0 [u]u [p\v 


K 

= 

0 

a[v]pu [u]pv a\p\uv 


L 


0 


( 40 ) 


for the jump functions defining u, u, and p. 

A nontrivial family of ODEs may be found by specifying the nonstan- 
dard Heaviside functions H', K' , and V if the determinant of the matrix of 
Equation (40) vanishes. The determinant of the matrix of Equation (40) is 
zero 


[u}u 

~[u]u 

0 


0 

[u]u 

w 

= 0. (41) 

a[u]pu 

[u\pv 

a \p]uu 



u 2 = 

7 pu. 

(42) 


Equations (41) and (42) imply that a nontrivial algebraic solution exists 
defining the ODEs specifying the nonstandard Heaviside functions for isen- 
tropic flow. Since shock propagation is not an isentropic process, Equation 
(40) will not be satisfied across a shock wave; therefore, Equation (40) is 
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replaced by the under- determined system 



l 

~ 

t 


r 

1 


H 


r -| 

[v}u 

— [u] v 0 


K 

= 

0 

0 

[u]u [p\v 




0 




L 




together with the entropy defined by 


( 43 ) 


s — c v ln(pz/ 7 ), 


(44) 


on either side of the shock wave, Zel’dovich and Raizer [23], p 53. 

The jump functions for u, u, and p are then determined by integrating 
matrix Equation (43) to determine the nonstandard Heaviside functions H, 
K, and L on an infinitesimal interval, (0, e), subject to boundary data 


H( 0) 


0 


m 


1 

K(0) 


0 

and 

K(e) 


1 

m 


0 


L{e) 


1 


(45) 


The boundary value problem defined by Equations (43) and (45) may 
be integrated to produce unambiguous relationships between the Heaviside 


25 



functions H, K, and L. To integrate H' and K', Equation (43) is rewritten 
in terms of a nonsingular matrix and an unknown function depending on L' 


u]u — [u\u 


H 

/ 

m 

0 [u]u 


K 


g(L') 


where 


m 


0 

g{L') 


-\p]vL' 


Equations (46) and (47) then imply 


~ 

/ 

- 

" 


- 

H 

1 

[u]u 

[u\u 


0 

K 

[v\ [u\v? 

0 

[u}u 


~W L ' 


(46) 


(47) 


(48) 


where 


H' 



(49) 


and 


K' 


W r t 
[u\u 


(50) 
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Then combining Equations (49) and (50) yields 


H' 



( 51 ) 


Integrating Equation (51) and applying the boundary conditions (45) shows 
that 


H « K, 


(52) 


with 


where 


= 1 

W\ui ’ 


Ui=Ui + c. 


(53) 


(54) 


Combining Equations (52) and (53) with Equation (50) produces 


K' = 



(55) 


And integrating Equation (55) with the boundary conditions (45) yields 


K w L, 


(56) 
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and 



( 57 ) 


Finally, by combining Equations (52) and (56) it follows that the microstruc- 
ture for the Heaviside functions for the specific volume, velocity, and pressure 
jump conditions are coincident across an inviscid shock wave 


H « _RT « L, (58) 

for an arbitrary infinitesimal interval (0,e). This result also holds if the 
equivalence relation « is replaced by = in Equations (45) and (58). The 
result of Equation (58) (with =) has been derived by Colombeau and Le 
Roux [10] using the theory of nonlinear generalized functions. 

Equation (58) was derived by manipulating products of the functions, H, 
K, and L and their derivatives. Because these functions satisfy (19), the 
mathematical structure applied in the analysis is not over-specified and can 
be used to derive the classical Rankine-Hugoniot jump conditions obtained 


by analyzing the problem in divergent form. For example, an expression for 



the shock speed c follows from Equations (53) and (54), which yields 

c= *? u iW\)- ( 59 ) 

Equation (59) is the shock speed obtained from the equations of motion in 
conservative form (which is the physically observed shock speed). A dis- 
cussion of the standard one- dimensional shock jump conditions is given in 
Chapter 1 of Gathers [25]. 

In case the relaxation phenomenon is important in the shock wave prop- 
agation process, the infinitesimal length scale over which the jump in the 
nonstandard Heaviside functions takes place should correspond to the phys- 
ical length over which of the slowest energy state reaches equilibrium. The 
physics of the relaxation phenomenon must be included in the differential 
equations to determine the microstructure of the Heaviside functions. 

3.3 Jump Conditions for Converging Shock Waves 

In analogy with the normal shock wave analysis of Section 3.2, the governing 
equations of motion are considered along characteristic curves in the space- 
time domain. Along the characteristic curves the equations of motion reduce 


29 



to ordinary differential equations. Following Zel’dovich and Raizer [23], the 
flow variables are assumed of the form 


p(0 = poR 2 ir(0,p(0 = pog(Oi and = R v (Oi ( 6 °) 


where the characteristic curves are defined by 

^ = R = A(—t) a ' 

In Equations (60) and (61), R (= R(t)), the location of the shock wave, and 
Po, the density in front of the shock wave, are the reference scales for the 
similarity solution; while A and a are constants. The functions of Equations 
(60) and (61) are defined on 



— oo < t < 0, R < r < oo, and 1 < £ < oo. (62) 


The converging shock fronts are located at £ = 1 and the shock speeds are 
given by D = R = aR/t. 

Recalling that shock wave propagation is not an isentropic process, Equa- 
tions (29) and (30) are used to define ODEs relating the density p, velocity 
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u, and pressure p across converging shock waves. Because the resulting dif- 
ferential equations will be under-determined, unique jump functions cannot 
be found using these equations. However, the analysis will yield unambigu- 
ous relationships between the jump functions for the flow parameters. And, 
as in the case of the inviscid normal shock, if one of the jump functions is 
specified, the other two jump functions can be determined. 

Along the characteristic curves, the similarity solutions n, g, and v are 
then assumed to be internal nonstandard jump functions across a shock front 


g(0 = go + b]G(£), 

(63) 

= Vo + [v]K(£), 

(64) 

tt( 0 = TTo + [7r]L(f), 

(65) 


In Equations (63), (64), and (65), each Heaviside predistribution is assumed 
to have its jump contained on the same interval (1, 1 + e), where e is an arbi- 
trary fixed infinitesimal. Here G, K, and L are assumed to be differentiable 
on (1, 1 + e), and G, K , L e *Li oc (R) as in Section 2.1. On the end points of 
the interval (1, 1 + e) the nonstandard Heaviside functions are subject to the 
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boundary data 


G(1) 


0 


G(l + e) 


1 

K( 1) 


0 

and 

K{ 1 + e) 


1 

Lil) 


0 


L(1 + e) 


1 


( 66 ) 


Behind a converging shock front, the density pi , velocity U \ , and pressure 
Pi, are assumed to take the limiting values for a strong shock 


Pi 


Po 


7 + 1 
7-1’ 


(67) 


ui 



( 68 ) 


and 

Pi = —T^PoD 2 , (69) 

7 + 1 

Thompson [24], p 495; where 7 is the ratio of specific heats. Combining the 
general functional form of the similarity solutions (60), with the nonstandard 
jump functions (63) to (65), the Heaviside boundary data (66), and the 
physical boundary data (67) to (69), yields the nonstandard jump function 
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parameters 


and 


5-0 = 1 and [5-] 


7 + 1 
7-1 


- 1 , 


no = 0 and [v 


2 

7 + 1’ 


(70) 

( 71 ) 


7r 0 — 0 and [it 


2 

7 + 1’ 


(72) 


Now to determine relationships between the nonstandard Heaviside func- 
tions, G, K and L, the equations of conservation of mass and momentum 
must be integrated. To this end, substituting Equations (60) into (29) and 
(30) produces the system of ODEs 


M 0 

0 [tt] 


K 

L 


f(G ) 

9(G) 


(73) 


where 


f(G ) 


f(hi 9)' - («/f + (ln#)> 

9(G) 


-0([(<* - l) / a]v + (v - £)v') 


(74) 


The system of Equations (73) have assumed that combinations of R, p and 
their time derivatives are functions of the constant a. Details of the deriva- 
tion of these equations and the scale factor identities are given in Chapter 
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XII of Zel’dovich and Razier [23]. 

The first Equation of (73) may be integrated to solve for either g or v. 
Since this equation forms a linear ODE in terms g, 


^ + (| + (lnj)> = ?(lnj)', 


(75) 


where (•)' = v is determined as a function g. Integrating Equation (74) 
on the interval (1, 1 + e) yields 




1 


C +1 9 (0- !-(« + !) / r K g(r)dr 


(76) 


To show that Equation (76) gives a jump function, it must be demon- 
strated that it satisfies the boundary conditions of Equations (66) and (71). 
From experimental shock data, g is assumed to be strictly monotonically 
increasing on the interval (1,1 + e); with this assumption the integral in 
Equation (76) satisfies 

0 < (« + 1) J T K g{r)d,T < {k + l)g(l + e) J r K dr. (77) 
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and 

W (1 + £)«-?-, (81) 
7 + 1 

so together with v(l) ~ 0. it follows from Equation (64) that v is a nonstan- 
dard jump function. Equation (80) is the homogeneous part of the solution 
of Equation (75), which implies that the nonhomogeneous term generates 
an infinitesimal contribution to the velocity jump function across the shock 
wave. 

Next, the second Equation in (73) gives a linear ODE for n in terms of v 
and g 

d'7T . CX 1 . . . . 

- = -g( v + (v-£) V . 82 

at, a 
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Integrating Equation (82) on the interval (1, 1 + e) produces 


vr(£) = ^ J ^ g(r)v(r)dT - J g (t)(v(t) - t)v'(t)cIt. (83) 

By combining Equations (76) and (83), 7 r is determined as a function of g. 
To simplify the resulting expression for n in terms of g, Equation (76) is 
replaced by Equation (80) in Equation (83), which yields 

^ dr + 

Using Equations (77) to (79), the first integral in Equation (84) is bounded 
by an infinitesimal, so that Equation (84) becomes 

7t(£) ~ J r-V(r)dr. (85) 

Integrating Equation (85) by parts then yields 

7 r(£) ~ ^ _K u(^) — u(l) + k J T~ K ~ 1 v(r)dT. (86) 

Note that the integral term in Equation (86) is bounded by an infinitesimal, 
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evaluating Equation (86) at £ = 1 and £ = 1 + e gives 


7r(l) ~ 0 and 7r(l + e) 


7 + 1 ’ 


(87) 


so 7r is a nonstandard jump function. And since Equation (78) implies that 


— « 1 
e ’ 


(88) 


on the interval (1, 1 + e), Equation (86) then becomes 


Jr(0 «»(£). 


(89) 


By combining Equations (63) to (65), (67) to (69), and Equation (89), it 
follows that the predistributions of the Heaviside functions G, K, and L for 
the density, velocity, and pressure jump conditions are in general distinct, 
across an inviscid shock wave 


Gy^K^L, 


(90) 
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and 


K « L, (91) 

for any infinitesimal interval (1, 1 + e). Moreover, the jump functions v, it 
and g are related by the integral equations 

v(() = ?W) + f^T-aMdr , (92) 

and 

J g(r)v(T)dT - j ^ g{r){v{r) - r)n'(r)dr. (93) 

The nonstandard functions g,v,ir G *Li oc (R) are solutions of the bound- 
ary value problem defined by Equations (73) and (74) and Equations (63) 
through (66) in the sense that the infinitesimal equivalence relation « re- 
places equality = in the equations of motion. 

4 Examples 

To exhibit examples for the microstructures of inviscid, converging shock 
waves, assume a density jump function g(£) e *Li oc (R) defined via Equations 
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(63) and (70) is given by 


9(0 = i + foK^h”, 


(94) 


on the interval (1, 1 + e), for n > 1. Also, assume the ratio of specific heats 
is given by 7 = 7/5 for all examples. 

Substituting Equation (94) into Equation (76) then yields 


u(0 = £ 


1 

£ K g(£) 


K + 1 

£ K g(0 


t k + t k \ 9 ]( 


T 



dr. 


(95) 


Fixing n — 1 and performing the integration in Equation (95) produces 


v(Z) = £ 


_j «+i ( {1 _ fck e +1 [ 9 ] e +2 

£ K g(£) £ K g(£) V e k + i e k + 2 

£ K g(0V eU + i + 


[g] 


c hi 2 


(96) 


By using Equations (66), (71), and (78) with Equation (96), v satisfies the 
boundary conditions and is a nonstandard jump function for velocity 


v(l) ~ 0 and v(l + e) « 


2 

7 + 1’ 


(97) 
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The jump function for pressure may be computed directly by combining 
Equations (94) and (96) with Equation (83) and integrating. However, the 
resulting form for the pressure contains many terms; so to simplify the anal- 
ysis, the pressure may be integrated numerically for the case of a converging 
spherical shock wave, k — 2. 

Figure 1 shows the density jump function specified by Equation (94) 
for n — 1. The shock wave jump interval (1, 1 + e) is normalized to the 
standard interval (0, 1) for all of the plots. The plots are generated with 
small real values for e and are used to convey the qualitative behavior of the 
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Figure 2: Computed velocity jump function for the density jump function 
shown in Figure 1: Circles e = 0.1, Boxes e = 0.01, Diamonds e = 0.001, 
Solid line limiting solution. 


nonstandard shock wave microstructure. 

Figure 2 shows the velocity jump function v computed from Equations 
(94) and (76). Approximations of the velocity are computed using Equa- 
tion (80) for the values of e = 0.1,0.01,0.001 and illustrate the asymptotic 
behavior the nonstandard jump functions. Moreover, as e approaches an in- 
finitesimal, the jump function for pressure approaches the jump function for 
velocity: n rj v, and hence, K r L. As a second example, the microstruc- 
ture is computed numerically for the case of a converging cylindrical shock 
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Velocity Jump Function Density Jump Function 




Figure 4: Computed velocity jump function for the density jump function 
shown in Figure 3: limiting solution e ~ 0. 
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wave, k — 1. From experimental data on normal shock waves, the field vari- 
ables have small variations near the upstream boundary of the shock layer, 
Muntz and Harnett [26]. To simulate this behavior, Equation (94) is con- 
sidered with n = 6, so that g\ 1) is small. Figure 3 shows the density jump 
function g on the normalized interval (0,1). The velocity jump function is 
then computed using Equation (76), and exhibited in Figure 4. Because the 
asymptotic behavior of the nonstandard Heaviside function for pressure im- 
plies that K ps L, it follows that Figure 4 also represents the pressure jump 
function 7r. 

The microstructure of converging shock waves depends on the geometrical 
and physical parameters, a, k, and 7. The similarity parameter a appears in 
an integral term of Equation (83) for the pressure jump function n. Because 
this term is bounded by an infinitesimal, a has little effect on the basic 
features of the inviscid shock microstructure and is neglected in the examples. 
Nonzero variation of the geometrical term k, which specifies the coordinate 
system (either 1 or 2), also has little effect on the gross microstructure. 
However, k was used to compute the examples and generates a small but 
noticeable variation for larger values of e. The ratio of specific heats 7 is the 
key physical parameter used to specify the macroscopic boundary conditions; 
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and as such, 7 fixes the gross features of the computed microstructure of 
converging shock waves for a given density jump function. 

5 Summary and Conclusions 

This paper derived shock wave jump conditions for one-dimensional, converg- 
ing shock waves in an inviscid, perfect gas using nonstandard analysis. Pre- 
distributions of the Heaviside functions were introduced and used to model 
the microstructure of the flow field across a shock wave governed by the clas- 
sical conservation equations. For flow variables with their jumps defined on 
the same infinitesimal interval, it was shown that the equations of motion 
written in nonconservative form produced unambiguous relations between 
the various nonstandard Heaviside functions. Because inviscid flow was as- 
sumed, the analysis only produced the relations between the field variable 
Heaviside functions. 

The present work was motivated by the authors’ interest in understanding 
how to add physical information to specify products of generalized functions 
contained in differential equations modeling physical problems with well de- 
fined regions of solution discontinuity. The main conclusions of the study are 
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as follows: 


1. Nonstandard analysis may be applied to identify predistributions of 
the Heaviside function defined on infinitesimal intervals that reproduce 
the standard jump function results from the theory of generalized func- 
tions. The predistributions considered here are piecewise differentiable 
elements of the function space *Li oc (R). 

2. Nonstandard asymptotic arguments were used to simplify the deriva- 
tion of the shock wave jump functions. 

3. The relationships between the internal Heaviside functions for one di- 
mensional, converging shock waves in an inviscid, perfect gas were de- 
rived unambiguously from the equations of motion in nonconservative 
form. Distinct nonstandard jump functions were obtained for the den- 
sity, velocity, and pressure in converging shock waves. 
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